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Abstract. We investigate the excitation spectrum of a model of N colour fermions 
with correlated hopping which can be solved by a nested Bethe ansatz. The gapless 
excitations of particle-hole type are calculated as well as the spin-wave like excitations 
which have a gap. Using general predictions of conformal field theory the long 
distance behaviour of some groundstate correlation functions are derived from a finite- 
size analysis of the gapless excitations. From the algebraic decay we show that 
for increasing particle density the correlation of so-called iV-multiplets of particles 
dominates over the density-density correlation. This indicates the presence of bound 
complexes of these iV-multiplets. This picture is also supported by the calculation of 
the effective mass of charge carriers. 
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1. Introduction 



The study of low-dimensional electronic systems with strong correlations has gained 
considerable importance due to the discovery of high-temperature superconductivity. 
The models which have been mostly studied are the one-dimensional Hubbard model 
[1] and the t — J model [2,3] at its supersymmetric point. These models are special 
because of their integrability [4-7] . A great deal of current research is directed to more 
extensive rigorous and exact results on correlated systems in one and two dimensions. 
However, in the latter case exact results are sparse. In contrast to this we have 
a different situation in one dimension where many models are integrable and show 
different physical behaviour. 

In this paper we continue the study of a fermionic model with correlated hopping 
[8]. In addition to on-site Coulomb interactions as in the above mentioned cases we 
are interested in modifications of the hopping terms of fermions in the vicinity of other 
particles. Unfortunately, such complex models are not integrable in general. However, 
for purely correlated hopping the situation looks much better. In [9] an integrable 
model of electrons was found which was generalized to fermions with arbitrary number 
N > 2 of colours [10]. The Hamiltonian of this model is given by 

L N / N \ 

H = - J2J2( c t,T C J+hr + c++i,t c j>) • exp -rj n J+e{T -T'),r' , (1) 

j = l t=1 V t' = i / 

where L is the length of the chain, c- T and c J;T are the creation and annihilation 
operators of fermions of colour r at site j (r = 1, N and j = 1, n J;T is the 

number operator and 9 is the step function 

T > t', 



lO, r<r'. 



(2) 



We employ the usual periodic boundary conditions cl+i, t = ci )T . Repulsion of 
particles corresponds to an interaction parameter rj > 0. In [8] the case iV = 2 of 
model (1) was investigated, however for negative values of rj. Note that rj and — rj are 
related by a particle-hole transformation. Using a Jordan- Wigner transformation the 
model can be formulated equivalently as a system of iV interacting XY chains [10]. 

In Sect. 2 we give a summary of the Bethe Ansatz equations and some groundstate 
properties. In Sect. 3 we study the excitations of the model. In Sect. 4 we derive the 



long-distance behaviour of correlation functions and compare the results with previous 
ones [8] for the special case of N = 2. In the Appendix a short-cut to the analysis of 
Bethe Ansatz equations is presented on the basis of inversion identities. 



2. The Bethe Ansatz 

The Bethe Ansatz for this model has been derived in [10] from which we quote the 
relevant equations. The eigenstates of the Hamiltonian are characterized by sets of 
wavenumbers kj, j = l,...,J\f, for M particles. There are additional Bethe Ansatz 
parameters (r = 1, N — 1 and a = 1, Mj^_ r ) the number of which is given 
by 

N-r 

M N _ r = E -A/}) (3) 
i=i 

where Afj denotes the number of particles of colour j. The parameters kj and A^ 
satisfy the set of nested Bethe Ansatz equations derived in [10] 



k 3 L+ E e (^- A i 1) ;|)=27r/ J , 

a=l 

E E (aw - a{T } ; I) - E ( A ^ r) - 4 } ; v) = ^ N ~ r) > 



with j = 1, A/"; r = 1, AT — 1; a = 1, M N _ r , and we have set A^ -* = kj. The 
phase shift function 6 is given by 



0(/c; 7]) = 2 arctan ^cothrytan ^/c^ , — 7r < < tt, 



(5) 



and Ij and Jq are integer (half- integer) numbers for odd (even) Mn-i + 1 and 
M r _i + J\f r+ i, respectively. Energy and momentum of the corresponding state are 
given by 

E = — 2 E cos 

(6) 

J = l \J = 1 f=l Q = l 



It follows by symmetry that the ground state corresponds to symmetric 
configurations with the same number of particles A/i for all colours i 

Mi = ^, 1 < z < N. (7) 

The ground state energy is calculated from (4)-(6) in the thermodynamic limit [10] 



-y- = -2 / cos kp(k)dk, (8) 

L J-K 

where the density function p(k) is determined as the solution of the integral equation 

,K 

2np(k) - cp(k- k')p{k')dk' = 1, 

J-K 

1 ~ smh\nr){N - 1)} , (9) 

^(fc) = 1 - — + 2 > exp(-n??) ±— — — -cosfn/c), 

v ; sinh(n?7iV) 

and X is determined by the subsidiary condition 



K M 

p(k)dk = — = p. (10) 

k L 



Equations (8)-(10) determine the ground state energy of the model as a function of 
the particle density p. The Appendix contains a derivation of eqs. (8)-(10) from 
(4)-(6) which is slightly different from the original treatment in [10]. 

In (7) we have assumed M a multiple of N . This in fact is the condition for a 
proper accomodation of the 'antiferromagnetic' ground state. For particle numbers M 
not multiples of N we always have a misfit resulting in 'spin-type' excitations which 
are treated below. This kind of excitations has a non-zero gap by which the bulk 
energies of the ground states of systems with Af not multiples of are increased. 



3. Excited states 

The system of Bethe Ansatz equations (4) admits many solutions depending on the 
choice of the parameters and Ja . In particular the solution for the ground state (8) 

(r) 

is characterized by a parameter set in which the Ij and Ja are consecutive integers 
(or half-integers) centered around the origin. It is expected that the low-lying states 
are obtained by small modifications to these configurations. 
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There are two types of elementary excitations leaving the particle number 
unchanged. The first one is of particle-hole type corresponding to raising one pseudo- 
particle in the ground state \kh\ < K to a higher level \k p \ > K. This excitation 
is obtained by a straight forward manipulation of the set of numbers Ij without 

(r) 

changing the parameters J& . This produces gapless excitations analogous to those in 
the antiferromagnetic spin-1/2 Heisenberg chain [11,12] and in the repulsive Hubbard 
model [13-15]. In the thermodynamic limit energy and momentum of these excitations 
are given by 



(11) 



p-po = Ep«)-EpW)» 

with energy-momentum dispersion of particle-hole excitations 

e(t?) = — 2 cos# + / sin kj(k)dk 

J-K 

r K 

p (&) = / j(k)dk. 

J-K 

The function j(k) is the solution of the integral equation 

rK 

2nj(k) - <p(k- k')j{k')dk' = $(Jfe - i?), 

J-K 



(12) 



(13) 



where $'(/c) = (p(k). Alternatively, the dressed energy function e can be obtained 
directly from an integral equation identical to the last one upon the replacement of 
j(k) and 3>(k — i?) by e(k) and eo(fc) = — 2cos/c, respectively. In general (13) can be 
solved only numerically. In certain limiting cases also analytic results can be obtained. 
For instance in the strong-coupling limit (rj — > oo) perturbation theory yields 



6(0) =-2 
p(0) 



1 + 



7V 



K - -sin2if 
2 



n 



TT-Kt 



+ 



2e~ 2 ^ sinK 
n- Kt 



(cos ■&+- (sin K - K cos A")) + O (e" 2?? ) , 
+ 0(e- 4 ^, 



. n tsmK n 
sin v H ?/ 

7T — A t 



N 



N 



(14) 

The excitations of the second type involve "r-holes" and "strings" in the sets 
A "r-hole" occurs when there is a jump in the sequence of Ja \ Physically it 



describes states where the colour of one bare particle has been changed. In a sense 
these excitations are of spin-wave type, but have a non-zero gap. The general energy- 
momentum excitation is given by N — 1 different dispersion functions 



JV-l vn- 



(15) 

p-po=y: E^ _r) H r) )- 



r=l j=l 

Each elementary excitation parametrized with a rapidity variable is given by 



r K 

e (N-r)^ = 2 / S [ n kj(k)dk, 

J ~ K (16) 



p (N-r)(#) = f j(k)dk, 
J-K 



where j(k) is the solution of an integral equation of the type (13), however with 
different right hand side 

r K 

2nj(k)- <p(k-k')j(k')dk' = ^ N - r \k-d), (17) 

J-K 



where 



{N - r){k) = + 2 f sinhMJV- r)] sin(nfc) 

N ^ smh(nr]N) n 

The lowest energies A( r ) of (16) are taken at t? = it. In Fig. 1 the gaps are shown 
for different values of iV and r\ in dependence of the particle density p. Not all 
combinations of the elementary excitations (16) are allowed, the numbers v r in (15) 
have to satisfy a selection rule 



N-l 



rv N - r = (mod N). (19) 



r=l 



For a derivation of (17) and (18) utilizing inversion identities see the Appendix. In 
general (16), (17) have to be solved numerically. In the strong-coupling limit (rj — > oo) 
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the analytic solution of (16)-(18) is 

2e- r7? / 1 



e (N-r) ^ = A (N-r) + ^_ - ± SU1 2K) (1 + COS 0) + O (e^ 2 ^) , 

(N-r)/q\ K N-r 2T lS mK / _ 4 \ 



A (JV-r) = 



,N — r tt-0 
iV 7T — Kt 
-2»? , 1 



^ _ r71 . n -2n JV —r-K — V . T ^ 

Ti = e r7? sin^-e 2?? — — ^sinK, ^0) 



yriV 



1 + ^— - i sin (l - T- V 2 "^) 
(sin K — K cos if) . 



In addition to simple holes excited states may also contain conjugate pairs of 
complex \ The additional "strings" do not contribute to energy and momentum 
which effect is very transparent in the inversion identity approach, see the Appendix. 
Such a cancellation phenomenon is well-known for the Heisenberg model, see for 
instance [12]. Concerning the spectrum of the model these states only lead to a 
degeneracy of the energy levels. With respect to the classification of all eigenstates, 
however, they do play an important role. 

Thus in the considered model there are no gapless excitations of spin wave type in 
contrast to the antiferromagnetic Heisenberg chain [11-12] and the repulsive Hubbard 
model [13-15]. Instead the analogous excitations correspond to collective motions of 
pseudo-particles which can be regarded as new particles with non-zero masses. 

The gap implies the binding of particles in multiplets of N particles with different 
colours. Any spatial separation of the constituents of these complexes would lead to 
two subsystems of the chain with particle numbers not multiples of N. As explained 
in Sect. 2 this results in an increase of the ground state energies by the order of 
the "spin-gap", and eventually leads to the exponential decay of the corresponding 
correlation function. 



4. Critical exponents of the correlation functions 

In the Bethe Ansatz approach it is a formidable task to deal with correlation functions. 
However, due to developments in two-dimensional conformal field theory the scaling 
dimensions describing the algebraic decay of correlation functions became accessible 
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[16,17]. According to this theory there is a one-to-one correspondence between the 
scaling dimensions and the spectrum of the quantum system on a (periodically closed) 
finite chain [18,19] 

2-kvf 



E j -E = Z -^( Xj+ N+ + N-), 

ii f \ (21) 

Pj -P = — (a, +N+ - N~J + 2dk F . 

Ej and Pj are the energy and momentum of the jth excited state, v F is the Fermi 
velocity. The Xj and Sj are the scaling dimensions and "spins" of the corresponding 
operators, iV + , N~ are non- negative integers and d is the number of particles excited 
from the left Fermi point to the right one. The ground state energy is expected to 
scale like [20,21] 

E = Le^ - ^fc, (22) 

where eoo is the ground state energy per site of the infinite system and c is the central 
charge characterizing the underlying conformal field theory. 

Consequently the scaling dimensions are obtained from the gaps due to finite-size 
effects in the spectrum of the Hamiltonian at criticality. In order to compute the 
finite-size corrections we use a method which has been developed by previous authors 
for the Heisenberg and Hubbard chains [22,23]. As a result we have obtained the 
ground state energy (22) with central charge c = 1 and the finite-size spectrum of the 
low-lying excitations 



E — En 



2 

AM 



+ [ti(K)] 2 d 2 + N + + N~ 



(23) 



2tt 

P = 2k F d+ — (dAAf + N+ - AT) , k F = irp. 



Here the dressed charge at the Fermi surface £,(K) is related in a simple way to the 
function p(k) in (9) 

Z(K) = 2np(K), (24) 

and AAf is the change of the particle number compared to the ground state. The 
non-negative integers N ± describe the excitations of particle-hole type in the vicinity 
of the Fermi points ±k F - Now we can read off all critical exponents from (23) and 
(21) 

(AM) 2 

x ^ d = ^ + m)]w + N+ + N-, (25) 

SAAf,d = dAM + N+ -N~, 
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where AJ\f and d can take the values 

AJ\f = 0, N, 2N, 

, AM ( 1 (26) 

generalizing the result in [24] for N = 2. The omitted integer values of AJ\f correspond 
to excitations with gap which are not related to correlations with algebraic decay. The 
two-point correlation functions of the scaling fields $a± t) with conformal weights 
A ± = (x± s)/2 are known to be [16,19] 

/ \ _ „ s . exp(— 2idkpx) 

(* A± (x, *)* i± (0, 0)) = { —^ { —L^ . (27) 

The dimensions of the descendant fields differ from A ± by the integers N ± . 

Let us first consider the density-density correlation function. In this case the 
correlation function is determined by excitations with unchanged number of particles 
AAf = 0. The leading contributions to the algebraic decay of the correlation are then 
given by d = 0, (N+ , N~) = (1,0) or (0, 1), and by d = 1/N, iV± = 0, respectively. 
This leads to the asymptotic form 

(p(r)p(O)) ~ pi + A ir ~ 2 + A 2 r~ a cos (Jjkpr^ , (28) 

where 

p(r)=f>, r , a (29) 

T = l 

Field-field correlators decay exponentially as the corresponding excitations change 
the particle number M by one unit, as remarked before. However, we may consider 
the correlation function of a iV-multiplet of fields, 

N 



0(r) = H Cr,r ■ (30) 



T = l 

In this case the particle number changes by AAf = N which according to (26) is the 
smallest non-zero value for which we have algebraic decay. We have d = or 1/2N 
for even or odd N, respectively. We thus obtain 

a f 1, f l/a, N even, 

O + r O ~ r ~P. J 13=1 (31) 

V 1 J 1 " \cos(*fr), M \ l/a + a/4, JV odd. 1 J 
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The characteristic exponent (29) is determined numerically and plotted in Fig. 2 for 
the cases AT = 3, 4. For N = 2 we refer to [8]. Analytically we find a(p = 0) = 2/N 2 
and a(p = N) = 2. This shows that for all A" there exists an interval \p c ,N] for 
which (3 < a. In this case the A^-multiplet correlation (31) has a slower decay than 
the density correlation (28) and thus dominates. This might indicate the presence of 
bound complexes of A^-multiplets in generalization of bound pairs for N = 2 which 
has been discussed in [8]. 

In order to get a more complete picture of the physical properties of the model 
we also investigate the conductivity and the effective transport mass. Following the 
ideas of [25,26] we study the Hamiltonian (1) with twisted boundary conditions with 
twisting angle (p. In this case the Bethe Ansatz equations read 

Mat- 



kjL+ e>(^- A i 1) ;|) = 2 ^ + ^ 



(32) 

E E e K r> - A Jr'; |) - E e ( A *' - = 2 -4 w -". 

(7 = ±1 0=1 /3=1 

The additional phase (p corresponds to an enclosed magnetic flux in the ring. 

The conductivity is directly proportional to the charge stiffness (Drude weight) 
D c which can be obtained from the change 

AE = D c ip 2 /L (33) 

in the ground state energy for small <p [25]. Using the charge stiffness we can define 

an effective transport mass m by the relation 

m D° c 

— = TT 34 

m e D c 

where D® = — sin (^) is the charge stiffness of the non-interacting (rj = 0) system 
and m e is the (bare) electron mass. 

Using the result (23) for the finite-size corrections with AAf = and d = cp/2Tt we 

find 

D c = ^v F e(K) = ±-J(K)Z(K). (35) 

Using the integral equation (9) and (24) we can calculate D c and m numerically for 
all densities p, see Fig. 3. For large particle densities p and strong interaction rj we 
observe an increase of the mass m by a factor A^ as compared to the non-interacting 
case. This is consistent with and supports our above picture that bound complexes of 
AT-multiplets are present in the system. 
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Appendix 



Here we describe the reduction of the Bethe Ansatz problem (4) involving kj and 
A^ parameters to equations involving only the kj variables. To this end we cite a set 
of equations which are equivalent to (4) and which appeared at an intermediate stage 
of the analysis in [10] 

e^ L = T 1 (L kj + |) . (36) 
T 1 (w) is the first of a set of functions T 1 (w), T 2 («),..., T N ~ 1 (u), satisfying the recursion 



T r («) = ' 



(37) 



q r (u 

with T N (u) = 1. The q- functions are given by 

q r (u) = J] sinh^-UWj, (38) 

in terms of the parameters A^ which are subject to the condition that all T r (u) be 
analytic functions with simple poles at ^Aa~^ — \. This condition is equivalent to 
the set of Bethe ansatz equations (4). Note that Aq = kj and M N = Af. 

On the right hand side of (37) one of the two summands usually dominates over 
the other. This provides the possibility to derive a simple functional equation for 
T 1 (it) in product form 

T\u)T\u + V)-T\u + (N - 1) V ) ~ 0( q °^'% vv (39) 

where the right hand side depends only on the kj parameters and the corrections to 
this relation are exponentially small in the thermodynamic limit. For the ground state 
(39) can be solved uniquely (apart from A^th roots of unity) under the condition of 
analyticity, absence of zeros in the physical regime — (N — l)rj <Re(w) < rj, and the 
pole structure at ^kj — \. The result can be written in the form 

M 

(40) 



where 



lnTl [\ k + 1) = %2 w J + 1 £ ® {k " h) ' 3 = °' 1; N ' 

^ ' 1 1=1 



w V 4^ n sinh(n?yA^) v y v y 



n=l 
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Using standard reasoning from this and (36) one derives the integral equation (9). 
Note that <p(k) = &'(k). The absolute ground state of the system corresponds to the 
choice of j = in (40). Other choices of j lead to excitations with energies of order 
1/L and give rise to a quantization of the d parameter in (26) in units of 1/N. 
For excitations of spin-wave type it is useful to define the excitation function 

T l (u) 

t(u) = lim — Yy (42) 

for which (39) directly implies 

t(u)t(u + rj)...t(u + {N — l)rj) = 1. (43) 

A first consequence of this functional equation is the double periodicity of the analytic 
continuation of t(u) 

t{u) = t(u + itt), 

(44) 

t( u ) = [t(u + rj)t(u + 2rj)...t(u + (N — 1)77)] 1 = t(u + Nrj). 

The elementary spin-wave excitations are characterized by pairs of zeros and poles 
which have to appear in pairs u z and u p . The distance of u z and u p has to be a 
multiple of 1] in order to lead to the cancellation (43) . We write 



u 



, P = l(v + i$Tr^), r = l,...,tf-l, (45) 



where t? is a rapidity variable and takes real numbers. The solution of (43) under the 
condition of analyticity and the zero/pole structure of (45) is 

lnt (^ + l) =^ (N ~ r) (k-$), (46) 

where the function on the right hand side was defined already in (18). Again by 
standard reasoning and (36) one finds the integral equation (17). It is noteworthy 
that in this approach the distribution of A-parameters is quite irrelevant. Different 
patterns, in particular simple holes with or without strings, lead to the same function 
(46). 
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Figure captions 

Figure 1. (a) Depiction of the gap of the massive excitations in dependence of the 
particle density p for N = 2 and different interaction strengths r\ = 0.5, 1, 2, oo. (b) 
The excitation gap for the four massive dispersions in the case N = 5 and r\ = 1. 

Figure 2. (a) Depiction of the critical exponent a for N = 3 and different values 
of r\ = 0, 0.1, 0.5 oo. For values of a above the broken line we have dominating 
iV-multiplet correlations, (b) The same for N = 4. 

Figure 3. The effective mass of the charge carriers for iV = 3 and different interaction 
strengths 7] = 0, 0.1, 0.5, 1, 10. For particle densities close to 3 and large interactions 
rj the mass approaches 3 times the bare mass m e . 
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